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Hadmn spectrum of QCD with one quark flavor 



1. Introduction and motivation 

QCD with one flavor of quarks {Nf = 1 QCD) radically differs from QCD with two or more 
flavors due to the absence of a chiral symmetry: the abelian symmetry of the one flavor theory 
is washed out at the quantum level by the Adler-Bell-Jackiw anomaly. Only a vector symmetry 
survives, related to the conservation of the quark number. As a consequence of this, the main 
features of the phase structure and mass spectrum of the single flavor theory strongly deviate from 
the familiar picture, affected by the spontaneous breaking of the non-abelian chiral symmetry, of 
ordinary QCD. These unphysical features of Nf = \ QCD explain the little attention reserved to 
this kind of setup in past simulations (however with the exceptions of ^ and 

This situation has changed in recent years, mainly due to the works of M. Creutz drawing the 
attention of the lattice community to open problems in the physical {Nf > 1) theory [Q, Since 
these aspects are not directly related to the spontaneous breaking of the chiral symmetry, they find 
an equivalent in the single flavor theory. The latter represents therefore a simple setup for their 
investigation. 

One question raised by Creutz |^, having a relevant phenomenological impact, is whether it 
is possible to define in an unambiguous way the case where one quark (say the u quark) becomes 
massless. The arguments against an unique definition of the massless limit [||] essentially rest 
upon the U{\) anomaly and should therefore hold a fortiori for the one flavor theory. A second 
aspect is the possibility of a spontaneous breaking of CP in QCD for special choices of the quark 
masses, conjectured for the first time by Dashen [^. According to the Vafa-Witten theorem [^] a 
prerequisite for the spontaneous breaking of a discrete symmetry is a non positive fermion measure, 
which in Nf = 1 QCD is possible for negative quark masses. The transition line is indeed expected 
to be located on the negative real quark mass axis in the extended complex parameter space. In 
the case of the multi-flavor theory, the transition is excluded for physical values of the quark mass, 
but its nearby presence can nevertheless affect numerical simulations on the lattice So the main 
features of this transition are not of academic interest only. 

Another intriguing aspect of one flavor QCD, emerging from string theory, is the connection 
with the supersymmetric Yang-Mills theory (SYM). The equivalence of the two theories in 

the bosonic sector [^] can be proven at the planar level of a particular large A^^ limit (orientifold 
large Nc limit) preserving balance between fermionic and bosonic degrees of freedom. Relics 
of SUSY are therefore expected m Nf = I QCD (with Nc = 3). A prediction of the orientifold 
equivalence [^], already studied in the literature [^, concerns in particular the size of the quark 
condensate. 

Another important place where relics of SUSY in Nf = 1 QCD can be investigated, considered 
in more detail in this contribution, is the low-lying bound-state spectrum [|T0|]. In SYM the mass 



patterns are strongly constrained by SUSY. In particular low-energy models [ ]11[ ] predict a low- 
lying chiral supermultiplet including two scalar particles with opposite parity. ' \n Nf = I QCD 
these two particles can be easily identified with the T] and the a meson (the former picking up 
a mass through the anomaly). On the basis of the planar equivalence, their mass ratio including 
0{\/Nc) corrections is expected to be rria/mri = Nc/ {Nc — 2) 



'For recent lattice simulations of ,yV=l SYM, see |T^, p3|]. 
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The Study of the mass spectrum of hadronic states requires reasonably large physical volumes, 
in order to be able to accommodate the bound-states in the finite box, and small (positive) quark 
masses. High statistics is required for a precise determination of the disconnected quark diagrams 
needed for the scalar meson masses, which are characterized by a high level of noise. High statis- 
tics is also important for the computation of the glueball masses. We apply here the Wilson lattice 
fermion action which has recently been shown [ [l5| , [l^ , |l7[ 18] to be well suited for such an in- 



vestigation. Preliminary results with Stout-smeared links [19] in the Wilson fermion action will 



be also presented. Following |18] we apply in the gauge sector the tree level improved Symanzik 
action (tlSym). 

The present exploratory study has been performed on 12^ • 24 and 16^-32 lattices with a lattice 
spacing corresponding in QCD units to a ~ 0.19fm and a ~ 0. 13fm, respectively. (We use the 



Sommer parameter [20] ro for setting the scale, fixed at the conventional value vq = 0.5 fm.) For 
the future we plan to run simulations closer to the continuum limit. 

As already mentioned, the sign of the quark determinant is an important issue mNf = \ QCD 
(in particular, a negative determinant triggers the CP-violating phase transition). In the continuum, 
the fermion determinant is positive for positive quaik mass. With Wilson lattice fermions for small 
quark masses, it can become negative due to quantum fluctuations. In most of our simulations the 
quark mass is large enough to prevent sign changes and the occurrence of a negative determinant is 
a rare event. For the lightest simulated quark masses however the sign of the quark determinant may 
potentially play a role and its impact in the hadron spectrum must be checked. In our simulations 
we could reach quite small quark masses down to lUq ~ 12MeV (m^ro ~ 0.03), corresponding to a 
pion mass mj^ ~ 270 MeV. 



As we have argued in [ IC ], it is useful to embed the Nf = 1 QCD theory in a partially quenched 
theory with additional quark flavors. A particularly symmetric choice consists in taking the {Ny) 
valence quark flavors degenerate with the sea quark: in this case the combined sea and valence 
sector is characterized by an exact S\]{Nv + 1) flavor symmetry. In this fictitious multi-flavor theory 
a PCAC quark mass can be naturally defined. We take this quantity as an operative definition of the 
quark mass for the (unitary) one-flavor theory. Also, a partially quenched chiral perturbation theory 
(PQChPT) can be set-up, exactly as in the Nf > 1 case. The latter reduces to an effective theory of 
the T] meson in the unitary sector without valence quarks. The predictions of this PQChPT will be 
compared against our numerical data. 

The plan of this contribution is as follows: in the next section the partially quenched viewpoint 
is introduced and PQChPT is considered for it. In Section ^ some information on the simulation 
algorithm and on the computation of the sign of the determinant are given. Section ^ is devoted 
to the presentation of our numerical results on the hadron spectrum, while Section || discusses the 
partially quenched data. The last section contains summary and outlook. 

2. Partially quenched QCD 

The symmetry of the one flavor theory can be artificially enhanced by adding extra valence 
quarks which are quenched, namely not taken into account in the Boltzmann-weight of the gauge 
configurations by their fermion determinants. A theoretical description of the resulting partially 



quenched theory can be obtained through the introduction of ghost quarks [21]. In this method the 
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functional integral over the ghost quark fields ij/ cancels the fermion determinant of the valence 
quarks i/Zy, 

det(7^D^+mv) 

and only the determinant of the sea (S) quark remains in the measure. In principle, one might 
consider any number of quenched valence quarks with any mass values. In our approach we take 
two valence quarks u and d with masses my and one sea quark s with mass m^. For our purpose 
the case of degenerate valence and sea quark mass my = nig is particularly convenient (which is 
admittedly an unconventional kind of partially quenching). Observe that in this symmetric setup 
the exact number of valence quarks Ny is immaterial, so our position NF=Nv+Nf = 2 + \ is 
just suggested by analogy with the case realized in nature. (Of course, in order to be able to build 
bound states containing two different quark flavors as mesons and nucleons, one needs Ny 

At the point of vanishing quark masses (see below) the generic partially quenched theory has 
a graded SU(A'^F|A^y)L S\]{Nf\Nv)r symmetry, which is broken spontaneously into a "flavor" 
symmetry SU(A^f |A^y), also vahd for non-vanishing degenerate quark masses. The SU(A^ir) sub- 
group represents the flavor symmetry in the combined sea and valence quark sectors. The latter 
symmetry implies that the hadronic bound states appear in exactly degenerate SU(A^ir) multiplets 
for my = nis- 

In particular, this extended theory contains a degenerate octet of pseudoscalar mesons ("pions" 
n°, a = 1, ... ,8) satisfying an SU(3)-symmetric PCAC relation. Considering the divergence of the 
axial-vector current A"^ and pseudoscalar density P" we can define the bare PCAC quark mass in 
lattice units as usual 

{dlA+p-) 

ampcAc = — ^ • (2.2) 

Here the indices + and — refer to the "charged" components corresponding to Xa ± iXb (with Xa,b 
some off-diagonal Gell-Mann matrices) and d* denotes the backward lattice derivative. Due to the 
exact SU(3)-symmetry, the renormalized quark mass corresponding to wjpcAc can be defined by an 
SU(3)-symmetric multiplicative renormalization: 

R _ ON 

As we will confirm numerically in sec. ^, the masses of the "pions" can be made to vanish 
by suitably tuning the bare quark mass on the lattice. In this situation the renormalized quark 



mass (|2.3[ ) vanishes, too. We stress here that the pions are not particles in the physical spectrum 
of the theory. Nevertheless their properties as mass and decay constant are well defined quantities 
which can be computed on the lattice. The same applies for the PCAC quark mass wip^^^, which can 
be therefore regarded as a potential candidate for a definition of the quark mass in this theory. 

2.1 Chiral perturbation theory 

The dependence of pion properties upon the quark masses can be determined in partially 



quenched chiral perturbation theory (PQChPT) [22, 23]. The effect of the finite lattice spacing 
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a can be also included [24, 25, 26, 27, 28]. The pseudo-Goldstone fields are parameterized by a 
graded matrix 

(2.4) 



U (x) = exp —^{x) 

which in our case is in the supergroup SU(3|2). (The normalization of Fq is chosen such that its 
phenomenological value is ~ 86MeV.) The commuting elements of the graded matrix represent 
the pseudo-Goldstone bosons made from a quark and an anti-quark with equal statistics, while the 
anticommuting elements represent pseudo-Goldstone fermions which are built from one fermionic 
quark and one bosonic quark. The supertrace of has to vanish, which can be implemented by a 
suitable choice of generators [29|. 

We have calculated both the masses and decay constants of the pseudo-Goldstone bosons in 
next-to-leading order of partially quenched chiral perturbation theory along the lines of Ref. [p^], 
including ff{a) lattice artifacts [^]. The quark masses enter the expressions in the combinations 



Xv = 2Bo m , Xs = 2Bo ms , ;i:pcac = 2Bo ml 



(2.5) 



with the usual leading order low -energy constant Bq ; the lattice spacing enters in the combination 

p = 2Woa, (2.6) 

where Wq is another, lattice-specific, low-energy constant. We have calculated the masses of the 
pions and mixed mesons (degenerate in the special case my = ms). The next-to-leading-order 
expression in terms of the (renormalized) PCAC quark mass is 



^PCAC + 



In 



(2L8-L5 + 2L6-L4);t;, 



2 

PCAC 



+ {Ws + W6-W5-W4- 2Lg +L5- 2Le + L4);t;pcAcP 



(2.7) 



where the usual next-to-leading order low-energy parameters L, appear, together with additional 
ones (Wj) describing lattice artifacts. For the decay constant we obtain in this case 



1 



Zpcac , Xpcac 



327r2f2 



In- 



A2 



+ [(L5 + U)Xpcac + {W5 + W4-L5- U)p] } . (2.8) 



Observe that as expected the results are independent of Ny- In particular, calculating the quan- 
tities in this model with Ny = 1, which corresponds to a representation the supergroup SU(2|1), 
reproduces (2/7) and (2.8). 

The analysis can be extended by relating the pion mass to the mass of the "physical" tj. The 
inclusion of the singlet can be achieved by relaxing the constraint of a vanishing supertrace [22, 2^], 
and associating it with the field 

^o{x) = sTr^{x). (2.9) 



The effective Lagrangian then contains additional terms depending on <I>o: 



(2.10) 
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where a and are free parameters in this context. We will use in the following the leading order 
expression for the mass of the T] , which reads 

2 ^% + ^PCAC 1 1 \ 

m, = ^^. (2.11) 
Our numerical results for allow to determine a and (see Section ^. 

3. Simulation 



For the SU(3) gauge sector we apply the tree-level improved Symanzik (tlSym) action [ |30| ] 
including planar rectangular (1x2) Wilson loops: 

S, = ^l,U t (l-^Ret/i^v^l+ci £ li-UeUlA] , (3.1) 

with ci = — 1/12 and normalization condition cq = 1 — 8ci. The fermionic part of the lattice action 
is the simple (unimproved) Wilson action. With the goal of improving the stability of the Monte 



Carlo evolution at small quark masses, we also started simulations with Stout-smeared links [ |190 
in the hopping matrix (see below). 



The update algorithm is a Polynomial Hybrid Monte Carlo algorithm (PHMC) [gjj, |32|] allow 



ing the simulation of an odd number of fermion species. The present version ||33|] is based on a 
two-step polynomial approximation of the inverse fermion matrix with stochastic correction in the 
update chain: a sequence of PHMC trajectories is followed by a Metropolis accept-reject step with 
a higher precision polynomial. The polynomial approximation scheme and the stochastic correc- 



tion in the update chain are taken over from the two-step multi -boson algorithm of Ref. p4|]. A 
correction factor C\U] in the measurement is associated with configurations for which eigenvalues 
of the (squared Hermitian) fermion matrix \U ] lie outside the validity interval of the polynomial 
approximation. We refer to [ ]T0| ] for more details on the algorithmic setup. 

As mentioned in the Introduction, the sign o\U] of the fermion determinant det2[[/] has also 
to be included in the reweighting of the configurations. The expectation value of a quantity A is 
therefore given by 

^ S[dU]a[U]C[U]A[U] 

^""^ ndU]a[U]C[U] ■ ^'-^^ 

For the computation of the sign o[U]we applied two methods. In the first we studied the spec- 
tral flow of the Hermitian fermion matrix [^]. For the K"-dependent computation of the low-lying 
eigenvalues of the Hermitian fermion matrix Q[U] we followed in this case Ref. [|36|]. Alternatively, 
we computed the (complex) spectrum of the non-Hermitian matrix concentrating on the lowest real 
eigenvalues: sign changes are signaled by negative real eigenvalues. We applied the ARPACK 



Arnoldi routines [ 37 ] on a transformed Dirac operator. The (polynomial) transformation was tuned 



such that the real eigenvalues were projected outside the ellipsoidal bulk containing the whole 



eigenvalue spectrum [p8[]. This allows for an efficient computation of the real eigenvalues [p^]. 
This latter method, on which we will rely in the future, delivers unambiguous results and can be 
simply automatized. 
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Table 1: Summary of the runs: 12^ • 24 and 16^ • 32 lattices have lowercase and uppercase labels, respec- 
tively. The bar indicates runs with Stout-link in the fermion action (see text). 







K 


^conf 


plaquette 


^nlaq 


ra/a 


a 


3.80 


0.1700 


5424 


0.546041(66) 


12.5 


2.66(4) 


b 


3.80 


0.1705 


3403 


0.546881(46) 


4.6 


2.67(5) 


c 


3.80 


0.1710 


2884 


0.547840(67) 


7.6 


2.69(5) 


A 


4.00 


0.1600 


1201 


0.581427(36) 


4.3 


3.56(5) 


B 


4.00 


0.1610 


1035 


0.582273(36) 


4.1 


3.61(5) 


C 


4.00 


0.1615 


1005 


0.582781(32) 


3.3 


3.73(5) 


A 


4.00 


0.1440 


5600 


0.577978(23) 


9.7 


3.74(3) 


B 


4.00 


0.1443 


5700 


0.578167(28) 


11.3 


3.83(5) 



3.1 Simulation details 

We performed simulations on a 12^-24 lattice with jS = 3.8 and on a 16^-32 with f5 =4.0. 
Information regarding the generated sets of configurations are reported in Table 

The sequences consisted of 3-6 PHMC individual trajectories. The precision of the first step 
of polynomial approximations was tuned such that the acceptance of the PHMC trajectories was 
about 0.80-0.85. The same acceptance was required for the Metropolis test by tuning the total 
length of the trajectory (1.5-1.8). This resulted in a relatively high total acceptance of 0.64-0.72. 
Optimization of the parameters of PHMC turned out to have a substantial impact on the integrated 
autocorrelation times of the average plaquette. 

In the case of a Stout-link we consider one step of isotropic smearing with p^y = p = 0.15, 
/I, V = 1, ... ,4. The Stout-smearing has in general the beneficial effect, compared to the unsmeared 
action, of reducing the fluctuations of the smallest eigenvalue of the (squared) hermitian matrix, 
with the result that less exceptional configurations are observed. This allowed us to obtain smooth 
simulations down to quite small pion masses niji ~ 270 MeV. 

Taking the values of r^/a at the highest K"'s for the runs at j8 = 3.8 and j3 = 4.0 and fixing 
ro = 0.5 fm by definition we obtain a = 0.186fm and a = 0. 134 fm, respectively. The extensions of 
the 12^ and 16^ lattices are roughly constant: L = 2.23 fm and L = 2.14fm. (The Stout-smearing 
leaves rg/a essentially unchanged.) 

For runs b, c, A and B there are cases where the eigenvalues of the fermion matrix are outside 
the approximation interval [s,X] and therefore C[?7] 7^ 1. In run c in particular there are 167 of such 
configuration out of 2884, 26 of them with negative sign. However, even in this case the average 
value a[?7]C[?7] is very near to one: 0.9842. The effect of the correction factors turns out to be 
quite weak in the case of the average plaquette and of ro/a: the effect on the average value of ro/a 
is only in the fifth digit (whereas the statistical error is in the third digit). This is not the case for 
low energy quantities as the low-lying hadron masses (see in the following). 
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Table 2: Results for light hadron masses xnNf = \ QCD (the result with the asterisk has been obtained with 
higher statistics: 4900). Note that the glueball masses were obtained at small time separations and hence 
could be overestimated (see also text). 









amQ++ 




a 


0.462(13) 


0.660(39) 


0.777(11) 


1.215(20) 


b 


0.403(11) 


0.629(29) 


0.685(10) 


1.116(38) 


c 


0.398(28) 


0.584(55) 


0.842(16) 


1.204(57) 


A 


0.455(17) 


0.607(57) 


1.083(79) 


1.006(15) 


B 


0.380(18) 


0.554(52) 


1.032(66) 


0.960(15) 


C 


0.316(22) 


0.613(67) 


0.943(41)* 


0.876(26) 



4. Hadron spectrum 

4.1 Mesons 

For the meson states we consider the simplest interpolating operators in the pseudoscalar and 
scalar sectors: 

77(0-): P{x) = y{x)y5¥{x) , (4.1) 
a(0+) : S{x) = i//(x)i/a(x) . (4.2) 

Corresponding states in the QCD spectrum are the 77'(958) and /o(600) (or a). In the case of the 
pseudoscalar mesons, invariance under the flavor group plays a special role when comparing with 
QCD states because of the U(l) axial anomaly. 

The disconnected diagrams of the hadron correlators of rj and a were computed by applying 
stochastic sources with complex Z2 noise and spin dilution. The method was already applied to the 
case of lattice SYM [p^]. In order to optimize the computational load, also considering autocor- 
relations, we analyzed typically every fifth configuration, with 20 stochastic estimates each. The 
resulting statistics is 400 — 600 on the smaller lattice and ~ 200 on the larger one. 

4.2 Baryons 

The simplest interpolating field in the baryon sector containing just one quark field is 

The low lying hadron state interpolated by the above operator is expected to have spin 3/2 and 
positive parity (|^). This corresponds to the A++(1232) of QCD if our dynamical fermion is 
interpreted as an u quark (the Q. baryon is more appropriate for larger quark masses). 



A difficulty arises since the Rarita-Schwinger spinor (4.3) also contains a spin 1 /2 component. 



We extract the wanted spin 3/2 component by projection pO|]: 

^3/2(0 = ^Tr [Gji{t)rjri + Gu{t)] , Gji{t) = £ {Aj{x,t)UO)) . (4.4) 
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''o'^PCAC 



Figure 1: The mass of the Ughtest physical particles in one-flavor QCD as a function of the bare PC AC 
quark mass. The masses are multiplied by the scale parameter ro in order to obtain dimensionless quantities. 
Open and full symbols refer to )3 = 3.8 and )3 = 4.0, respectively. 



Since the baryon correlator does not contain disconnected diagrams, our full statistics could be 
taken for the computation of the masses in this case, namely 3000^000 on the smaller lattice and 
~ 1000 on the larger one. 



4.3 Glueballs 

Spin states are also projected by purely gluonic operators. These are the glueballs, a well 
known object of investigation in lattice QCD. In particular the 0++ glueball has the same quantum 
numbers as the a meson. In this first investigation we neglect possible mixings between the two 
states and consider only diagonal correlators. 

We used the single spatial plaquette to obtain the mass of the 0++ ground state. To increase the 
overlap of the operator with this state we performed APE smearing and also applied variational 
methods [ p^ ] to obtain optimal glueball operators from linear combinations of the basic operators. 



4.4 Results 

The results for the hadron masses (only available for the runs without Stout-smearing) are 
reported in lattice units in Table In Fig. |l| the hadron masses are plotted as a function of the 



bare PCAC quark mass /WpcAc ( |2.3D defined in the partially quenched picture. Since we use physical 
units here, results from the two lattice spacings can be compared. The scaling is satisfactory for 
the case of Tj , whose mass could be computed with the best accuracy. The determination of the a 
meson mass seems to require large statistics. 

The effect of the sign of the determinant in the hadron spectrum was investigated by computing 
the masses with or without the inclusion of the sign factor in the reweighting procedure. Only in 
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Table 3: The PCAC quark mass mpcAc, the pion mass niji and non-renormalized decay constant f^, and the 
nucleon mass in lattice units (only for the runs without Stout-smearing). 







aruji 


afn 




a 


0.02771(45) 


0.3908(24) 


0.1838(11) 


1.0439(54) 


b 


0.01951(39) 


0.3292(25) 


0.1730(15) 


0.956(27) 


c 


0.0108(12) 


0.253(10) 


0.156(10) 


1.011(51) 


A 


0.04290(36) 


0.4132(21) 


0.1449(9) 


0.9018(44) 


B 


0.02561(31) 


0.3199(22) 


0.1289(10) 


0.7978(53) 


C 


0.01700(30) 


0.2635(24) 


0.1188(12) 


0.734(10) 


A 


0.01532(34) 


0.2316(49) 


0.09747(15) 




B 


0.00886(75) 


0.1994(74) 


0.0852(49) 





the case of our run at the lightest quark mass, run c, a sizeable effect can be observed: here the sign 
of the determinant pushes up the masses by 7 — 10%. 

We observe that our statistics is not large enough to obtain an accurate estimate of the glue- 
ball masses. In particular, the results reported in Table |2| could be overestimated. Indeed, due to 
the high level of noise, large time-separations could not be included in the determinations; it is 
therefore possible that the latter are contaminated by excited states. In order to enhance the statis- 
tics we decided to store the gauge configuration more frequently (as was already applied for the 
continuation of run C). 



5. Partially quenched analysis 

The results for the partially quenched sector are collected in Table || and shown in Figure 0. 
This also includes the nucleon mass (only for the runs without Stout-smearing). 

The partially quenched ChPT formulae are used to extract the corresponding low-energy co- 
efficients from the pion data. Considering the number of lattice data at our disposal, a full fit 
including all the terms in the ChPT formulae is not possible, so we take only the continuum terms 
into account. We fitted the data for both j3 values simultaneously neglecting the dependence of the 
renormalizations factors Za and Zp upon the lattice coupling constant. Introducing the one-flavor 
low-energy constants 

A3 = 47rfoexp{647r2(L4 + L5-2L6-2L8)}, 

A4 = 47rFoexp{647r2(L4 + L5)}, (5.1) 
the fit formulae of the renormalized values reduce to 

™2 _ I Zpcac ] Zpcac 
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® 



<$> 



N, p=3.8 
71, p=3.8 
N, p=4.0 
71, p=4.0 



@ 



0.05 0.1 

''o'^PCAC 



0.15 



Figure 2: The mass of the valence pion and nucleon as a function of the bare PCAC quark mass. Open and 
full symbols refer to j3 = 3.8 and j3 =4.0, respectively. 



fn _ 1 XvCKC ^PCAC 



FoV2 327r2Fo2 a2 

The data and the fitted curves are shown in Fig. ^. 

In Oder to improve the numerical results for the universal low-energy constants A3 4, which do 
not explicitly depend on the lattice spacing a, we also performed fits to the ratios [43, 15] 



2 ' f 



(5.3) 



For this calculation we restricted ourself to the data at jS =4.0 with reference point atK" = 0.1615. 
We obtain in this case the results 

— = 10.0 ±2.6, (5.4) 
^0 

— = 31.5 ±14.3, (5.5) 
^0 



which, interestingly, are compatible with phenomenological values obtained from ordinary QCD [fl4[]. 
The errors are however quite large (we hope to improve these determinations in the future). 

In addition, we investigated the relation between the mass of the pion and of the physical f] , 
reducing to formula ( 2.1 1 ) at leading-order. For this purpose we fitted simultaneously m\ and 
as a function of the PCAC quark mass, again considering only j3 = 4.0. This yields to 

a = -0.03(19), am<i. = 0.18(8), (5.6) 

suggesting a vanishing a. Using the value of ro/a extrapolated to vanishing PCAC quark mass and 
setting a = we find 

amc = 0.19(2) or romo = 0.72(10) , (5.7) 
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Figure 3: Pion masses squared and pion decay constants in lattice units and the results of the PQChPT fit. 



which means 

m<,> = 284±40MeV (5.8) 

in physical units. 



The value of mo can also be obtained from the Witten-Veneziano formula [ 45 1 



ml = j^^Xr (5.9) 



valid at leading-order in the ('t Hooft) large 'Nc limit. An estimate of the quenched topological 
susceptibility present in the hterature is Xt = (193 ib9MeV)^ [|4^. Using our value for which 
is subject to a sizeable statistical error, one would obtain = 450 ± 170 MeV. 

6. Summary and outlook 

This first Monte Carlo investigation of Nf = I QCD reveals the qualitative features of the low 
lying hadron spectrum of this theory. The lightest hadron is the pseudoscalar rj meson (see Table ^ 
and Figure |I|) while the scalar meson, the a, is about a factor 1.5 heavier. It is interesting to compare 
our data with the estimate in [ p7| ] nifj/m-q ^ Nc/ {Nc — 2) = 3 for Nc = 3. The above prediction 
applies for the massless theory and one could expect the agreement to improve for smaller quark 
masses. Our bare quark masses (estimated from the PCAC quark mass in the valence analysis) 
range between 10 MeV and 60 MeV, while the lightest pion mass is ~ 270 MeV. 

The lightest baryon, the A (|^), is by about a factor 3 heavier than the rj meson. The lightest 
scalar mass obtained with a glueball 0++ operator lies between the a meson and the A baryon 
mass. However, this mass could be overestimated, since, due the high level of noise, only small 
time-separations could be included in the analysis. 

In general, the mass measurements have relatively large errors between 3-10%. In order to 
obtain more quantitative results, larger statistics and smaller quark masses are required. We hope 



to be able to make progresses in both directions |48] with our new simulations using Stout-smeared 
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links in the fermion action. Some preliminary results were already presented in this contribution 
(see [49] for a test of this formulation in twisted mass QCD with Nf = 2). 

The introduction of a partially quenched extension of the single flavor theory with valence 
quarks allows to define the bare quark mass in terms of the PCAC quark mass of the fictitious 
multi-flavor theory. The computation of the bare quark mass is intricate in the unitary theory due 
to the absence of a chiral symmetry (the arguments of [^] regard the definition of a renormalized 
quark mass). Comparison of lattice data with partially quenched chiral perturbation theory allowed 
the determination of some of the low-energy constants of the chiral Lagrangian. The latter are 
compatible, even if with large error, with recent lattice determinations for Nf = 2 QCD. 

A further direction of investigation for the future [18] is the CP- violating phase transition 
expected at negative quark masses [^. For this aspect of the single flavor theory the non-positivity 
of the fermion measure plays an essential role. 
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